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Abstract 

The behavior of fermions with two spin states that interact with 
a large scattering length is constrained by universal relations that 
hold for any state of the system. These relations involve a cen- 
tral property of the system called the contact, which measures 
the number of pairs of fermions with different spins that have 
small separations. The contact controls the thermodynamics of 
the system as well as the large-momentum and high-frequency 
tails of correlation functions. This review summarizes the current 
theoretical and experimental status of these universal relations. 



to be published as a chapter in BCS-BEC Crossover and the 
Unitary Fermi Gas (Lecture Notes in Physics), edited by Wil- 
helm Zwerger (Springer, 2011). 
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1 Introduction 

Particles with short-range interactions that produce a large scattering length 
have universal properties that depend only on the scattering length pQ. A 
system consisting of such particles is strongly interacting in the sense that 
there are effects of the interactions that must be treated nonperturbatively. 
These strong interactions give rise to strong correlations among the particless. 
Many theoretical methods, even if they are nonperturbative, are inadequate 
for dealing with such strong correlations. However, such a system is also 
governed by universal relations that follow from the short-distance and short- 
time dynamics associated with the large scattering length. These universal 
relations provide powerful constraints on the behavior of the system. They 
hold for any state of the system: few-body or many-body, ground state or 
nonzero temperature, homogeneous or in a trapping potential, normal state 
or superfluid, balanced in the two spin states or imbalanced. They connect 
various properties of the system, ranging from thermodynamic variables to 
large-momentum and high-frequency tails of correlation functions. 

The systems for which the universal relations have been most extensively 
studied are those consisting of fermions with two spin states. The universal 
relations that have been derived thus far all involve a property of the system 
called the contact, which measures the number of pairs of fermions in the 
two spin states with small separations. Many of these relation were first 
derived by Shina Tan, and they are known as the Tan relations [21 [31 E] - 
Tan derived these relations by exploiting the fact that the large scattering 
length can be taken into account through boundary conditions on the many- 
body Schrodinger wavefunction for otherwise noninteracting particles. The 
universal relations can also be derived concisely within a quantum field theory 
framework [S] , where they follow from renormalization and from the operator 
product expansion. Such a framework facilitates the derivation of additional 
universal relations and the systematic inclusion of corrections associated with 
the nonzero range of the interactions. 

In this review, we summarize the current theoretical and experimental 
status of universal relations for systems consisting of fermions with two spin 
states and a large scattering length. We begin in Section [2] by presenting 
the Tan relations. In Section [3| we discuss the physical interpretation of the 
contact and we provide some illustrative examples. In Section [4j we present 
other universal relations that have been derived more recently. In Section [5j 
we describe exciting recent developments in the field of ultracold atoms in- 
volving measurements of the contact and experimental tests of the universal 
relations. In Section [6j we discuss the derivation of the universal relations, 
with an emphasis on quantum field theory methods. 
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The Tan relations were derived by Shina Tan in a series of three papers 
[3j H]. The first two of these papers were written in 2005, but they were not 
published until 2008, when all three papers were published back-to-back in 
Annals of Physics. 

The Tan relations apply to systems consisting of fermions with two spin 
states whose scattering length a is large compared to the range tq of their 
interactions. We will refer to the fermions as atoms and label the two spin 
states by an index a with values 1 and 2. In a many-body system, the number 
densities n a and the temperature T must also be small enough that the 

— 1/3 

corresponding length scales are large compared to the range: n a ^> tq 
and At 3> ?"q, where At = {2nh 2 /niksT) 1 / 2 . If the system is in an external 
trapping potential V(r) = \muj 2 r 2 , the length scale associated with the trap 
should also be large compared to the range: ih/mw) 1 / 2 ^> Tq. 

The Tan relations involve an extensive quantity, the contact C, which is 
the integral over space of a local quantity, the contact density C{R): 

C= Jd 3 RC{R). (1) 

We proceed to present the Tan relations in chronological order. 



2. 1 Tails of distributions 

In the first of Tan's two 2005 papers, he derived three universal relations [2]. 
The first was for the tails of the momentum distributions n a {k) for the two 
spin states a = 1,2: 

Tail of the momentum distribution. The distributions of the wavevec- 
tor k have power-law tails at large k: 

n a {k) — > C/k 4 . (2) 

The coefficient C is the contact and it is the same for both spin states. 

The asymptotic behavior in Eq. |2]) actually applies only in the scaling region 
| cz | 1 <C k <C r$ . The wavenumber must also be larger than the scales asso- 
ciated with the system, such as n 1 / 3 , A^ 1 , and (muj/K) 1 / 2 . The momentum 
distributions in Eq. ^ have been normalized so that the total number of 
atoms in the spin state a is 



* - / (2vr )3 n *( fc )- 



(3) 
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The universal relation in Eq. ([2]) implies that the contact is positive defi- 
nite and has dimension (length) _1 . Thus the contact density has dimensions 
(length)- 4 . 

The total energy E of the system is the sum of the kinetic energy T, the 
interaction energy U, and the energy V associated with an external potential: 

E = T + U + V. (4) 

The kinetic energy T (which should not be confused with the temperature) 
can be expressed as an integral over the momentum distribution: 




The asymptotic behavior of the momentum distribution in Eq. ([2j) implies 
that T is ultraviolet divergent. This divergence actually occurs only in the 
zero-range limit ro — > 0. For interactions with a finite range, the integral in 
Eq. ([5l is cut off by the range and therefore has a contribution that behaves 
like l/r as r — > 0. Thus the physical interpretation of the ultraviolet diver- 
gence is that T is sensitive to the range. The second Tan relation in Ref. [2] 
implies that the sensitivity of the kinetic energy to the range is cancelled by 
the interaction energy: 

Energy relation. The sum of the kinetic and interaction energies is ul- 
traviolet finite and it is completely determined by the momentum distri- 
butions n a (k) and the contact C: 

x ^ r d 3 k h 2 k 2 ( .,, c\ h 2 „ 

In the integral on the right side, the subtraction term cancels the tail of the 
momentum distribution and makes the integral convergent in the ultraviolet. 
The sum of the two terms in Eq. ^ proportional to the contact C is the 
interaction energy. The first of those two terms is ultraviolet divergent. Thus 
the interaction energy is sensitive to the range, but that sensitivity is exactly 
cancelled by the kinetic energy. The last term in Eq. ^ is the interaction 
energy that remains after subtracting the divergent term. Remarkably, it is 
also proportional to the contact. 

The third Tan relation in Ref. [5] gives the asymptotic behavior of the 
correlation function for the densities of the two spin states at short distances: 

Density-density correlator at short distances. The correlation be- 
tween the number densities for the two spin states at points separated by 
a small distance r diverges as 1/r 2 and the coefficient of the divergence is 
proportional to the contact density: 
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( ni (R+ §r) n 2 (R- |r)> JL (A - A) C (H). (7) 

Tan also pointed out that the contact density appears in the short-distance 
expansion for the correlator of the quantum field operators that create and 



annihilate the atoms. This expansion will be discussed in Section 6.4.4 



2.2 Changes in the scattering length 

From the three universal relations described above, one might conclude that 
the contact is an esoteric property of the system that has only to do with tails 
of distributions. In the second of Tan's 2005 papers [3J, he derived another 
universal relation that makes it clear that the contact is an absolutely central 
property of the system: 

Adiabatic relation. The rate of change of the energy due to a small 
change in the inverse scattering length is proportional to the contact: 



dE \ h 2 



da 1 J a 4nm 



a (8) 



The derivative is evaluated with the entropy S held fixed. The particle 
numbers N\ and N% are also implicitly held fixed. 

In the simplest case, E is just an energy eigenvalue. The adiabatic relation 
also holds for any statistical mixture of eigenstates if the derivative is eval- 
uated with the occupation numbers held fixed. By the adiabatic theorem 
of quantum mechanics, if the scattering length changes sufficiently slowly 
with time, the occupation numbers remain constant. Thus if the contact C is 
known as a function of a, the adiabatic relation in Eq. ^ can be integrated 
to obtain the accumulated change in E. 

The adiabatic relation can also be expressed in terms of the derivative of 
the free energy F = E — TS with the temperature T held fixed: 



dF \ Ti 2 



da 1 / rp Aum 



a (9) 



As pointed out by Tan, this implies that the contact determines the ther- 
modynamics of the system. Given the contact of a system as a function of 
the scattering length a and other variables, such as N\, N 2 , and T, the free 
energy F can be obtained by integrating Eq. ^ with respect to a. A conve- 
nient boundary condition is provided by the limit a — ¥ 0~ , in which the atoms 
are noninteracting. From F, one can determine all the other thermodynamic 
functions. 
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If one uses Eq. ([2]) to define the contact in terms of the tail of the mo- 
mentum distribution, this appears to be a case of the tail wagging the dog. 
The thermodynamic behavior of the system seems to be determined by the 
tail of the momentum distribution. However the proper interpretation is that 
the contact is a central property of the system that determines both the 
thermodynamics and the tail of the momentum distribution. 

The adiabatic relation in Eq. {[sj) determines the change in the total en- 
ergy when the scattering length changes very slowly. Tan also considered the 
opposite limit in which the scattering length changes very rapidly [3]: 

Sudden change in the scattering length. If the scattering length is 
changed suddenly from a to a', the change in the total energy is propor- 
tional to the contact: 

*E—£-(±-±)c, do) 
47rm \a' a J 

where C is the initial value of the contact. 

This result requires the time scale for the sudden change in scattering length 
to be much slower than the time scale mr^/h associated with the range. Tan 
also presented a more general result for the change in the energy due to a 
rapid change in the scattering length, which will be described in Section |4.5[ 



2.3 Additional Tan relations 

In Tan's 2008 paper, he derived two additional universal relations that apply 
for specific forms of the external potential [3]. 

Virial theorem. For a system in a harmonic trapping potential, the com- 
ponents of the energy E in Eq. Q satisfy 

H 2 

T + U -V = -- C. (11) 

oTTma 

The virial theorem for the unitary limit a — ±oo was first derived and also 
verified experimentally by Thomas, Kinast, and Turlapov [BJ. The virial the- 



orem in Eq. (11) is the generalization to finite scattering length. 



Pressure relation. For a homogeneous system, the pressure and the en- 
ergy density are related by 

V= 2 -£ + -^C. (12) 
6 iZirma 
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The pressure relation was actually first derived in Ref. [3] for the special case 
of a balanced gas in which the two spin states have equal populations. The 
derivation was extended to the general case in Ref. [3J. 

If there are inelastic 2-body scattering processes with a large energy re- 
lease, they will result in a decrease in the number of low-energy atoms. Tan 
realized that the rate at which the number density of low-energy atoms 
decreases is proportional to the contact density C [7]. The proportionality 
constant was first given in Ref. 5 . If there are inelastic 2-body scattering 
channels, the scattering length a has a negative imaginary part. The propor- 
tionality constant in the universal relation can be expressed in terms of that 
complex scattering length: 

Inelastic 2-body losses. If there are inelastic 2-body scattering pro- 
cesses with a large energy release, the number density of low-energy atoms 
decreases at a rate that is proportional to the contact density: 

d h(— Ima) „._, 

Jt n ^ R) = -^W c{R) - (13) 



3 What is the Contact? 

Given the Tan relations described in Section [2] it is evident that the contact 
is a central property of the system. But what is it? In this section, we provide 
an intuitive interpretation of the contact. We also provide additional insights 
into the contact by giving analytic expressions in some simple cases. 



3.1 Intuitive interpretation 

An intuitive interpretation of the contact density can be derived from the 
universal relation for the density-density correlator in Eq. ([7]). That relation 
can be expressed in the form 

(m(R + ri )n 2 (R + r 2 )> — > 2| 1 |9 C(fl). (14) 

107T \ri — r 2 \ 

If we integrate both t*i and r 2 over a ball of radius s, we obtain 

N pail (R,s)^^C(R). (15) 

The left side simply counts the number of pairs inside that ball, which is the 
product N1N2 of the number of atoms in the two spin states. The volume of 
that ball is V = f^s 3 - One might naively expect the number of pairs to scale 
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as V 2 as V — > 0. However, according to Eq. (15), it scales instead as V 4 / 3 . 
That scaling behavior applies only for s smaller than \a\ and also smaller 
than the length scales associated with the system, such as n" 1 / 3 , At, and 
(%/mu)) 1 / 2 . The scaling behavior extends down to s of order the range tq. 

A naive definition of the density of pairs is the limit as V — > of N\ N2/V 2 , 
where N± and N2 are the numbers of atoms in the volume V. This quantity 
has dimensions (length) ~ 6 . The result in Eq. (15) implies that the combina- 
tion with a nontrivial small-volume limit is JV1JV2/V 4 ' 3 , which has dimensions 
(length) . Thus a more appropriate definition of the local pair density is the 
small-volume limit of N1N2/V 4 / 3 , up to a normalization constant that can 
be chosen by convention. The unusual dimensions of (length) ~ 4 for this lo- 
cal pair density can be expressed concisely by saying that this quantity has 
scaling dimension 4. The difference —2 between the scaling dimension and 
the naive dimension 6 is called the anomalous dimension. This anomalous 
dimension comes from the strong correlations associated with the large scat- 
tering length. This anomalous scaling behavior implies that the number of 
pairs in a very small volume is much larger than one would naively expect by 
extrapolating from larger volumes. The contact density is a measure of the 
local pair density that takes into account this anomalous scaling behavior. 

Further intuition for the contact can be gleaned from the universal relation 
for the tail of the momentum distribution in Eq. ([2]). It implies that the 
number of atoms in either spin state whose wavenumber k is larger than K 
is 

> K) = ^C, (16) 

provided K is in the scaling region l/\a\ <C K <C 1/ro and is larger than the 
wavenumber scales set by the system. Thus the contact is a measure of the 
number of atoms with large momentum. 

When interpreting the contact density as a measure of the local pair den- 
sity, one should not confuse those pairs with Cooper pairs, which are pairs 
with a specific momentum correlation. Under conditions in which Cooper 
pairs are well defined, the typical separation of the atoms in a Cooper pair 
is much larger than the interparticle spacing. The number of Cooper pairs in 
a volume V that is comparable to or smaller than the volume per particle is 
not well defined. It certainly does not have the anomalous scaling behavior 
V~ 4 / 3 of the total number of pairs. 



3.2 Few-body systems 

The adiabatic relation in Eq. ([8| can be used as an operational definition of 
the contact. If the energy of a system is known as a function of the scattering 
length, we can simply differentiate to get an expression for the contact. A 
simple example in the case a > is the weakly-bound diatomic molecule, 
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or shallow dimer, consisting of two atoms with spins 1 and 2. The universal 
result for its binding energy is h 2 /ma 2 . Its energy is therefore 



h 2 



(17) 



Using the adiabatic relation in Eq. ([8]) , we find that the contact for the dimer 
is 

Cdimer = ■ (18) 

a 

Blume and Daily have calculated the contact numerically for the ground 
state of 4 fermions, 2 in each spin state, trapped in a harmonic potential 
and interacting through a potential with a large adjustable scattering length 
a [HJ. They determined the contact as a function of a using four different 
universal relations: the tail of the momentum distribution in Eq. the 
density-density correlator at short distances in Eq. |7]) , the adiabatic relation 



in Eq. (|8|), and the virial theorem in Eq. (111. The small differences between 
the four determinations of the contact were compatible with systematic errors 
associated with the nonzero range of the potential. 



3.3 Balanced homogeneous gas 

The contact density C for the homogeneous gas can be obtained by using the 
adiabatic relation in Eq. ([8]) as an operational definition. Dividing both sides 
of the equation by the volume, the relation can be expressed in the form 



c= Wd£ 
fi 2 da 



Analytic results for the energy density £ are available in various limits, and 
they can be used to obtain analytic expressions for the contact density. 

We first consider the case of a balanced gas, in which the two spin states 
are equally populated, at zero temperature. The total number density n = 
2ni — 2n 2 determines the Fermi momentum: kp = (S^n) 1 ^ 3 . The ground 
state is determined by the dimcnsionless interaction variable 1/kpa. As this 
variable changes from — oo to to +oo, the ground state changes smoothly 
from a mixture of two weakly-interacting Fermi gases to a Bose-Einstein 
condensate of diatomic molecules. The ground state is always a superfluid. 
The mechanism for superfluidity evolves smoothly from the BCS mechanism, 
which is the Cooper pairing of atoms in the two spin states, to the BEC 
mechanism, which is the Bose-Einstein condensation of dimers. 

In the BCS limit a — >• 0~, the energy density can be expanded in powers 
of kpa: 
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h 2 k% ( 10 



Using Eq. (19 1, we find that the contact density in the BCS limit is 

C — > Air 2 n 2 a 2 . (21) 

This is proportional to the low-energy cross section Ana 2 and to the product 
{n/2) 2 of the number densities of the two spin states. The contact density 
decreases to as a — > 0~ , but it decreases only as a 2 . This result emphasizes 
that the contact density is not closely related to the density of Cooper pairs, 
which decreases exponentially as exp(— 7r/fci?|a|) as a — > _ . 

In the unitary limit a —¥ ±oo, the energy density can be expanded in 
powers of l/kpa: 

c h 2 k% (, c 



107r 2 m V kpa 



i , ••• • (22) 



where £ and £ are numerical constants. Using Eq. (191, we find that the 
contact density in the unitary limit is 

C — ► ^(37r 2 n) 4 / 3 . (23) 

Since the interaction provides no length scale in the unitary limit, the contact 
density must be proportional to kp by dimensional analysis. An estimate 



( 1 for the numerical constant in Eq. ( 23 1 can be obtained from numerical 
calculations of the energy density near the unitary limit using quantum Monte 
Carlo methods [TU] . A more precise value can be obtained from numerical 
calculations of the density-density correlator in the unitary limit using the 
fixed-node diffusion Monte Carlo method [TT] . Using the Tan relation for the 
density-density correlator in Eq. ([7]), we obtain the value 

C « 0.95. (24) 

In the BEG limit a — > + , the energy density can be expanded in powers 
of (kpa) 3 / 2 : 



^ _ h 2 n TTh 2 n 2 a dd 
2ma 2 Am 



where a d d ~ 0.60 a is the dimer-dimer scattering length [T5| . The leading 
term is the total binding energy density for dimers with number density n/2 



and binding energy given by Eq. (17). The second term is the mean-field 



energy of a Bose-Einstein condensate of dimers with dimer-dimer scattering 



length a d d- Using Eq. (19), we find that the contact density in the BEC limit 
is 

Ann 

C — > . (26 

a 
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Fig. 1 The dimensionless contact density s = C/kp for the balanced homogeneous gas as 
a function of the dimensionless coupling strength 1/kpa, from Ref. 1131 . The left dashed 
line is the leading contribution in the BCS limit given by Eq. | |21| . The right dashed line 
is the leading contribution in the BEC limit given by Eq. j26j. The contact density in the 
unitary limit may be underpredicted by about 16%. 



This is equal to the contact 8n/a for a dimer, which is given in Eq. (181, 
multiplied by the dimer number density n/2. The contact density diverges 
as 1/a as a — > + . The first correction to the leading term in Eq. ( 26 1 is 
suppressed by a factor of (kpa) 3 . 

The contact density C for the ground state of the balanced homogeneous 
gas can be expressed as k F multiplied by a monotonically-increasing dimen- 
sionless function of l/kpa. Haussmann, Punk, and Zwerger have calculated 
this function numerically [13] . They used the Luttinger-Ward self-consistent 
formalism to calculate the single particle spectral functions. The contact den- 
sity was determined using the Tan relation in Eq. ^ for the tail of the mo- 
mentum distribution. Their result for the contact density as a function of 
1 / kpa is shown in Fig. [I] It interpolates smoothly between the BCS limit in 
Eq. (21 1 and the BEC limit in Eq. (26). Their result in the unitary limit cor- 
responds to a value ( ss 0.80 for the constant in Eq. (23). This is about 16% 
smaller than the value in Eq. ( 24 ) that was obtained from a fixed- node diffu- 
sion Monte Carlo calculation. This difference reflects a systematic theoretical 
error in the calculational method of Ref. [T3] ■ 

The dependence of the contact density on the temperature has been de- 
termined analytically in various limits by Yu, Bruun, and Baym [14] . In the 
low-temperature limit, the leading thermal contribution to the contact den- 
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sity comes from phonons. The leading thermal correction at small T increases 
like T 4 [TJJ. It becomes significant when T is comparable to the Fermi tem- 
perature: ksTp = % k F /2m. In the BCS limit and in the unitary limit, the 
thermal contribution differs from the contact at T = by a factor of (T/Tp) 4 
multiplied by a numerical constant. In the BEC limit, the thermal contribu- 
tion is suppressed by a factor of (T/Tp) 4 and by a further factor of (kpa) 1 ^ 2 . 

The contact density at high temperature can be calculated using the virial 
expansion. The high-temperature limitQwas determined by Yu, Bruun, and 
Baym [13] : 

8n 2 h 2 n 2 

Since the contact density increases as T 4 at small T and decreases as 1/T at 
sufficiently large T, it must achieve a maximum somewhere in between. The 
maximum occurs for T of order Tp . The maximum is pronounced only when 
a is near the unitary limit [T3] . 

Palestini, Perali, Pieri, and Strinati have calculated the contact density for 
the balanced homogeneous gas at nonzero temperature using a diagrammatic 
t-matrix approximation [TS]. They determined the contact from the high- 
momentum tail of the momentum distribution in Eq. ([2| and from the high- 
frequency tail of the radio-frequency transition rate (which is discussed in 
Section 

Thus far, the calculational methods that have been used to calculate the 
contact density numerically involve uncontrolled approximations. While they 
may be accurate in certain limits, there may also be regions of kpa and T/Tp 
in which the systematic theoretical errors are not negligible. 



3-4 Strongly -imbalanced homogeneous gas 

We now consider the strongly-imbalanced gas, in which a tiny population of 
minority atoms in state 2 is immersed in a system of atoms in state 1. The 
minority atoms can be considered as a dilute gas of impurities in the Fermi 
sea of majority atoms. In the homogeneous gas with number densities ri\ and 
ri2, the ground state is determined by the dimensionless interaction variable 
1/kpia, where kp\ = (Q^ni) 1 ^ is the Fermi wavenumber for the majority 
atoms. In the BCS limit a — > 0~, the impurity particle is an atom in state 2. 
In the BEC limit a —¥ + , the impurity particle is the dimer whose binding 



energy is given by Eq. (17). Using a diagrammatic Monte Carlo method, 
Prokof'ev and Svistunov have shown that there is a phase transition at a 
critical value a c of the scattering length given by 1/fc.F id c = 0.90 ± 0.02 [16 . 
As 1/kpia increases through this critical value, the impurity changes from a 

1 In Ref. 1141 . the contact density was denoted by 16ir 2 h 2 C and the total number density 
was denoted by In instead of n. 
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Fig. 2 The dimensionless contact density s = C/(67r 2 fcj? 1 ri2) for the strongly-imbalanced 
homogeneous gas as a function of the dimensionless coupling strength 1/kFia, from 
Ref. |17| . The dashed line is the leading contribution in the BEC limit given by Eq. \29\ . 
The dotted and dash-dotted lines are continuations of the solid line past the phase transi- 
tion into metastable regions. The position of the phase transition may be underpredicted 
by about 7%. 



quasiparticle associated with the atom in state 2, which is called a polaron, 
to a dimer quasiparticle. 

Analytic expressions for the contact density for the ground state of the 
strongly-imbalanced homogeneous gas can be obtained from the energy den- 



sity using Eq. (19). In the BCS limit a — > , the contact density is 

C — > 167r 2 nin 2 a 2 . (28) 

This is proportional to the low-energy cross section 47ra 2 and to the product 
ri\ni of the two number densities. The contact density is an increasing func- 
tion of 1 1 a that is smooth except for a discontinuity at the phase transition 
at a — a c . In the BEC limit a — > + , the contact density is 

C 87rn2 . (29) 



This is equal to the contact 8ir/a for the dimer, which is given in Eq. (18), 
multiplied by the dimer number density n 2 . 

The contact density C for the ground state of the strongly-imbalanced 
homogeneous gas can be expressed as kp\ni multiplied by a monotonically- 
increasing dimensionless function of 1/fcpia. This function has been calcu- 
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lated numerically by Punk, Dumitrescu, and Zwerger |17j . They used the 
adiabatic relation in Eq. ( 19 1 to obtain the contact from the ground state 
energy density, which they calculated using a variational method that gives a 
fairly good approximation to the results from the diagrammatic Monte Carlo 
method [H]. The results of Ref. !!?] for the contact density as a function of 
1/fcpia are shown in Fig. [2j For large negative values of l/kp^a, the contact 
density can be approximated by the BCS limit in Eq. (28). It is predicted to 
increase to about 5 kp\n 2 in the unitary limit a — > ±oo and then to about 
20 kp\ni at the phase transition a c , where it is predicted to have a disconti- 
nuity of about 20%. The prediction of the variational method for the position 
of the phase transition is (/cfiOc) 1 ~ 0.84, which is about 7% smaller than 
the diagrammatic Monte Carlo result 0.90 ± 0.02. This difference reflects a 
systematic theoretical error in the variational method of Ref. [T3] • For large 
positive values of 1/kpia, the contact density can be approximated by the 
BEC limit in Eq. f29l. 



4 Other Universal Relations 

Many new universal relations involving the contact have been discovered in 
recent years. They reveal that the contact plays a central role in many of the 
most important probes for ultracold atoms. 



4-1 rf spectroscopy 

Given a system of atoms in spin states 1 and 2, a radio-frequency (rf) signal 
that is tuned to near the transition frequency between an atom in state 2 
and an atom in a third spin state 3 can transform the atoms in state 2 
coherently into linear combinations of atoms in states 2 and 3. These atoms 
can be subsequently transformed by decohering processes into a mixture of 
atoms in states 2 and 3. The net effect is a transfer of atoms from the state 
2 to the state 3. The transition rate r(uj) for this process depends on the 
frequency ui of the rf signal. It is convenient to choose the offset for uj to be the 
transition frequency for a single atom. The transition rate for an extremely 
dilute sample of N 2 atoms is then a delta function at oj = : 

r(u) — > 7ri? 2 8(u) N 2 , (30) 

where fl is the Rabi frequency associated with the transition. In a many- 
body system consisting of atoms in states 1 and 2, r(ui) can be modified by 
initial-state interactions between atoms in states 1 and 2 and by final-state 
interactions between atoms in state 3 and atoms in states 1 or 2. However 
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the effects of these interactions are constrained by a sum rule [HI Qjj] : 

/oo 
duj r(u) = ttL? 2 N 2 . (31) 
-co 

If the atoms interact through large pair scattering lengths = a, 013, and 
023, there are universal relations that govern the rf transition rate r(ui). One 
of these universal relations is a sum rule derived by Punk and Zwerger 1201 
and by Baym, Pethick, Yu, and Zwierlein [21] : 

r°° %q 2 / 1 1 \ 

/ d WW r( W )=— ) Ci 2 , (32) 

J -00 4m \a 12 a 13 J 

where Gxi = C is the contact for atoms in states 1 and 2. The term propor- 
tional to l/ai3 comes from final-state interactions between atoms in states 1 



and 3. If we divide the sum rule in Eq. (32) by the sum rule in Eq. (31 1, we 
get an expression for the frequency shift (uj) averaged over the system. This 
frequency shift is called the clock shift. The universal relation for the clock 
shift has several interesting features. If 013 = a 12 , the clock shift vanishes 
because of a symmetry relating atoms 2 and 3. The clock shift has smooth 
behavior in the unitary limit ayi — > ±00. This behavior was first observed in 
experiments on rf spectroscopy in 6 Li atoms |22j . If we take the limit 013 — > 



in Eq. (32), the sum rule diverges. This implies that if the scattering length 
a 13 is not large, the clock shift is sensitive to the range. 

Another universal relation for rf spectroscopy is that the high-frequency 
tail of r(u>) is proportional to the contact. The general result for large scat- 
tering lengths ai2 and 013 was derived in Ref. 23J: 

r(w) -> a ( ^W275 2 T 7^ Ci2 - (33) 

4:iru:(mui/n) 1 ' (a 13 +mui/n) 



The asymptotic behavior in Eq. ( 33 ) holds when w is much larger than the 
many-body frequency scales hkp/m and hsT/h, but still much smaller than 
the frequency scale U/mrQ associated with the range. If lo ^> fi/ma\ 3 , the 
high-frequency tail decreases as w~ 5 / 2 . The result if the scattering length 013 
is not large can be obtained by taking the limit a\ 3 —¥ in Eq. ( |33[ ): 

In this case, the high-frequency tail decreases as uj^ 3 ^ 2 . This scaling behavior 
was derived in Ref. [23] and the coefficient was first calculated correctly by 
Schneider and Randeria |25j . The scaling behavior was also pointed ou t in 



Ref. [35]. If r(u>) decreases asymptotically as uj 3 / 2 , the sum rule in Eq. (32 ) 



diverges. Thus it is the high-frequency tail in Eq. (34) that makes this sum 



rule sensitive to the range in the case of 013 that is not large. 
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4-2 Photoassociation 



Photoassociation uses a laser to transfer pairs of low-energy atoms into an ex- 
cited molecular state with very high energy. The wavefunction of the molecule 
has support only over very short distances much smaller than the range To, 
so the pair of atoms must be very close together to have a reasonable prob- 
ability of making the transition. If there is a closed-channel molecule near 
the 2-atom threshold that can be excited by the laser, it can dominate the 
photoassociation rate. The rate is then proportional to the number N mo i of 
closed-channel molecules: 

Q 2 

r = — N mol , (35) 
7 

where fl is the Rabi frequency of the laser and 7 is the line width of the 
excited molecule. 

Werner, Tarruel, and Castin [28] and Zhang and Leggett [29] pointed out 
that if the large scattering length a comes from the tuning of the magnetic 
field B to near a Feshbach resonance associated with this closed-channel 
molecule, then N mo \ is proportional to the contact. If the zero of energy is 
chosen to coincide with the threshold for the atoms, the rate of change of the 
energy of the system with respect to the magnetic field can be expressed as 

7S = -Mmoi N moh (36) 
dis 

where /i mo i is the difference between the magnetic moment of the closed- 
channel molecule and twice the magnetic moment of an atom in the open 
channel. The scattering length a(B) near a Feshbach resonance at Bq can be 
parametrized by 

a(B) = a bg (l - ■ (37) 



Combining Eqs. (36) and (37) and using the adiabatic relation in Eq. ([8]), 
we obtain an expression for the number of closed-channel molecules that is 
proportional to the contact 28J: 

R A 2 

Nmo1 = 4nlA-(B-B )] 2C > (38) 

where i?* is a positive length that characterizes the width of the Feshbach 
resonance: 

h 2 

R* = ~r • (39) 

m[x mol a hg A 

That length can also be expressed as i?* = —\r s , where r s is the effective 
range at the center of the resonance B — B . 



Universal Relations for Fermions with Large Scattering Length 



IV 



The universal relation between the number of closed-channel molecules 
and the contact was previously derived formally by Braaten, Kang and Plat- 
ter [3D], but they did not make the connection to photoassociation. 



4-3 Structure factors 

Structure factors encode information about density-density correlations in a 
system. For a many-body system of fermions with two spin states, the correla- 
tions between the densities of the two spin states are particularly important. 
The corresponding static structure factor Si2(?) for a homogeneous system 
is the Fourier transform in the separation vector ri — r 2 of the correlator 
(ni(ri)n 2 (r 2 )) of the two densities. The dynamic structure factor Si2(q,uj) 
is the Fourier transform in the separation vector and the time interval be- 
tween the two densities. It encodes information about the degrees of freedom 
that can be excited by density fluctuations. The static structure factor can 
be obtained by integrating S , 12 (<7,w) over oj. 

If the scattering length is large, the static structure factor Si 2 (q) has a 
high- momentum tail that decreases like 1/q [51] . Hu, Liu, and Drummond 
have pointed out that this tail is proportional to the contact density [52] : 

SMq)^U---^)c. (40) 
8 V q iraq^ J 



The normalization of Si 2 (g) m Ref. Z2\ differs from that in Eq. (40) by a 
factor of 2/N, where N = Ni+N 2 is the total number of atoms. The universal 
relation in Eq. ( |40[ ) follows simply by taking the Fourier transform of the Tan 
relation in Eq. Wj for the density-density correlation at short distances. 

Son and Thompson have studied the dynamic structure factor S±2(q,u) 
in the unitary limit [33] . They showed that the leading contribution in the 
scaling limit oj — } oo and q — > oo with x = hq 2 /2muj fixed is proportional 
to the contact density. The coefficient of C is (mw 3 /^)" 1 / 2 multiplied by 
a complicated function of the dimensionless scaling variable x, which they 
calculated analytically. For small x, their result reduces to 

S 12 (q,w)-+ 2 V C. (41) 

Taylor and Randeria have also determined the high-frequency tail of the 
dynamic structure factor [34) . Their result for the limits q — > followed by 



oj — > oo differs from that in Eq. (41 ) by a factor of 3/2. 
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4-4 Viscosity spectral functions 



Taylor and Randeria have derived universal relations for the viscosity spectral 
functions of a homogeneous gas [34] . They found that the shear viscosity 
spectral function has a high-frequency tail that is proportional to the 
contact density: 

^ -> m^&jm c (42) 

They also derived a sum rule for r](uj): 

h 2 C \ irh ( „ 3h 2 



duJ y 1 ^ - ifewp J = t r " To^ C J ' (43) 

where £ is the energy density. Given the high-frequency tail of r][pj) in 



Eq. (42), the subtraction term on the left side of Eq. (43) is necessary to 
make the integral convergent. Enss, Haussmann, and Zwerger obtained a re- 
sult for the high-frequency tail of rj(ui) that also decreases as cj" 1 / 2 but has 



a coefficient smaller than that in Eq. ( 42 ) by a factor of 2 /3 [35] . 



Taylor and Randeria have also derived a sum rule for the bulk viscosity 
spectral function ((uj) [51] : 



?i 3 fdC 



where the derivative is taken with the entropy per particle held fixed. Since 



the spectral function £(u>) is positive definite, the sum rule in Eq. (44 ) implies 
that the contact density C is a strictly increasing function of a -1 . 



4-5 Rapid change in the scattering length 

In Ref. [3J, Tan showed that if the scattering length a(t) is time dependent, 
the total energy of the system changes at a rate that is proportional to the 
contact: 

^E(t) = -£ n C(t)b { t), (45) 

where b(t) = l/a(t) and C(i) is the instantaneous contact at time t. If the 
external potential V(r) is also changing with tim e, t here is an additional 



term proportional to V(r) on the right side of Eq. (45). Tan referred to that 
equation as the dynamic sweep theorem. 

The simplest case of a time-dependent scattering length is a sudden change 



in a, for which the change in the energy is given in Eq. ( 10 1 . Tan also presented 



a more general result for a scattering length a(t) that changes rapidly enough 
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that the contact does not have time to evolve significantly from its original 
value. If the scattering length a(t) changes over a short time interval T from 
an initial value a(0) to a final value a(T), the change in the total energy is 

AE=-^—(^--^- + J^- fdt fdt'Vt^b{t)b{t') \ C, 
Anm \ a(T) a(0) V mn J Q ' ' 



(46) 



where C is the contact at the initial time t = 0. The time interval T must 
be short compared to the time scales for the evolution of the system. It must 
also be long compared to the time scale mr^jfi set by the range. 

Son and Thompson have also considered rapid changes in the scattering 
length for a system that is initially in the unitary limit |33j . The simplest 
case is a small-amplitude oscillation of the inverse scattering length about 
the unitary limit: a(i) _1 = 7ocos(wt). If to is large compared to the Fermi 
frequency hkp/2m 7 the energy density increases at a rate that is proportional 
to the contact density: 

dt 8nm V m ^ ^ 

The frequency cu must also be small compared to the frequency scale h/mr^ 
set by the range. 

Another case of a rapidly varying scattering length considered in Ref. [3"5] 
is a pulse a(£) -1 that differs from zero only during a short time interval. The 
change in the total energy is proportional to the contact: 

*5/2 / poo \ 

^=^37*01 do, VE? |7(«)| a ) C, (48) 
where 7(0;) is the Fourier transform of the pulse: 7(0;) = J^dt exp(iw<:)a(<) _1 . 



This is consistent with the more general result in Eq. (46), which does not 
require the unitary limit before and after the pulse. 



5 Making Contact with Experiment 

There are some exciting recent developments in the study of the universal 
relations using experiments with ultracold atoms. They involve various mea- 
surements of the contact and the experimental verification of some of the 
universal relations. 
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5.1 Photoassociation 



The Hulet group at Rice University in Houston has measured the photoassoci- 
ation rate of a balanced mixture of 6 Li atoms in the lowest two hyperfine spin 
states [35] • The scattering length was controlled by using the Feshbach reso- 
nance near 834 G. They measured the photoassociation rate at various values 
of the scattering length, with 1/kpa ranging from about —1.4 to about +5.4. 
The photoassociation laser can excite the closed-channel molecule that is re- 
sponsible for the Feshbach resonance. The photoassociation rate can therefore 
be interpreted as a measurement of the number of closed-channel molecules. 
According to the universal relation in Eq. (38), this is proportional to the 
contact. However, at the time of the experiment, the concept of the contact 
was still unfamiliar. The number of closed-channel molecules was expected to 
be proportional to the square of the order parameter |zi| 2 [36]. In the BEC 
limit, |/i| 2 oc hp/a, which has the same dependence on a as the contact in 
Eq. (26). In the BCS limit, the order parameter decreases exponentially with 
l/|a|: |Z\| 2 oc fcp exp(— w/(kir\a\)). This behavior is dramatically different 
from the contact in Eq. (21), which decreases like a 2 as a — >• 0~. Neverthe- 
less, the measurements of the photoassociation rate were compatible with 
the assumption that it was proportional to \A\ 2 , even for negative values of 
1/kFd as large as —1.4. 

The first analysis of the data from Ref. in terms of the contact was 
carried out by Werner, Tarruell, and Castin [28 . Their results are shown in 
Fig. [3] The contact extracted from the measured number of closed-channel 
molecules was in reasonable agreement with a theoretical prediction using the 
local density approximation, with the contact density for the homogeneous 
system obtained by interpolating between the BCS, unitary, and BEC limits. 



5.2 Static structure factor 

The static and dynamic structure factors for systems consisting of ultracold 
atoms can be probed by using Bragg spectroscopy. Bragg scattering is a two- 
photon process in which an atom absorbs a photon from one laser beam and 
emits a photon into a second laser beam. The net effect is the transfer of a 
selected momentum Tik and a selected energy Huj to the atom, where k is the 
difference between the wavevectors of the two lasers and ui is the difference 
between their frequencies. 

The Vale group at Swinburne University of Technology in Melbourne has 
used Bragg spectroscopy to study the static structure factor 5i2(fc) [37] . 
They used a balanced mixture of 6 Li atoms in the lowest two hyperfine spin 
states. The scattering length was controlled by using the Feshbach resonance 
near 834 G. In Ref. [35], they reported measurements of the static structure 
factor Si2(k) as a function of a for k = 4.8 kp. These measurements are 
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l/(k p a) 

Fig. 3 The derivative of the energy E with respect to 1/a determined from measurements 
of the photoassociation rate of a trapped gas of 6 Li atoms, from Ref. 1281 . The data come 
from measurements of the number of closed-channel molecules in Ref. 1361 . The solid line is 
a theoretical prediction using the local density approximation, with the contact density for 
the homogeneous system obtained by interpolating between the BCS, unitary, and BEC 
limits. The symbol + indicates the prediction for the contact in the unitary limit. The 
dashed lines are extrapolations from the BCS and BEC limits. 



in good agreement with the universal relation for the large-momentum tail 
in Eq. (40), with the contact density for the homogeneous sytem obtained 
by interpolating between the BCS, unitary, and BEC limits. In Ref. [38] . 
they also reported measurements of S12 (fc) as a function of k for l/kpa — 
—0.2, 0, and +0.3, which are shown in Fig. [4j The measurements are linear 
in kp/k, as predicted by the universal relation in Eq. (40). The slope is 
predicted to be proportional to the contact C. The contact for their trapped 
system at zero temperature was calculated in Ref. [38] using the local density 
approximation and a below-threshold Gaussian fluctuation theory for the 
homogeneous system. For l/kpa = —0.2 and +0.3, the slope agrees well 
the universal relation. In the unitary limit \/kpa = 0, the slope is smaller 
than predicted. The discrepancy could be attributed to the effects of nonzero 
temperature. 
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Fig. 4 The structure factor S(k) for a trapped gas of 6 Li atoms as a function of kp/k 
for three values of 1/kpa, from Ref. | 38| . The solid lines are the universal predictions from 
Eq. ( |40| using the contact obtained from a theoretical calculation. The dashed line is a 
linear fit to the data for 1/kpa = 0. 




5.3 Comparing measurements of the contact 



The Jin group at JILA in Boulder has measured the contact C for a trapped 
gas of atoms using three independent methods [55]. They used a balanced 
mixture of 40 K atoms in the two lowest hyperfine spin states at a temperature 
of about 0.1 Tp. The scattering length was controlled by using the Feshbach 
resonance near 201 G. 

The first method for measuring C used the Tan relation in Eq. ([2| for the 
high- momentum tail of the momentum distribution n 2 (fc). The interactions 
between the trapped atoms were turned off by changing the magnetic field to 
209 G where the scattering length vanishes. The trapping potential was then 
turned off, and the momentum distribution was measured from the ballistic 
expansion of the cloud of atoms. The contact C is the large-momentum limit 
of k 4 n2(k). It was measured for values of 1/kpa ranging from about —2.7 to 
about +0.2. 

The second method for measuring the contact used the universal relation 



in Eq. (34) for the high-frequency tail of the radio-frequency (rf) transition 



rate. The rf signal was used to transfer atoms from state 2 into a third spin 
state 3 for which the pair scattering length 013 is not large. The rf transition 
rate -T(w) was determined by measuring the number of atoms transferred. 
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(k F a)' 1 

Fig. 5 Three measurements of the dimensionless contact C/kp for a trapped gas of 40 K 
atoms as a function of 1/kpa, from Ref. |39l . The Fermi wavenumber kp is defined by the 
Fermi energy for the trapped system: Ep = h 2 kp/2m. Two of the data sets are from the 
tail of the momentum distribution measured directly by ballistic expansion (solid dots) 
and indirectly by photocmission spectrometry (open dots). The third data set is from the 
high-frequency tail of the rf lineshape (stars). 

The contact was then determined from the behavior of r(u) at large u. It 
was measured for values of 1/kpa ranging from about —1.7 to about +0.2. 

The third method for measuring C also used the Tan relation in Eq. ([2]), 
but the tail of the momentum distribution was determined by photoemission 
spectroscopy (PES). This involves using momentum-resolved rf spectroscopy 
to measure the distribution n2{k,ui) of the momentum and energy of atoms 
in state 2, and then integrating over to to determine the momentum distri- 
bution n 2 {k). The contact C is the large- momentum limit of k i n2(k). It was 
measured for the same values of 1/kpa as the second method. 

The three sets of measurements of the contact by the Jin group are 
shown in Figure |] The results fr om the three methods are all consistent. 
They lie close to the theoretical prediction of Ref. j28j , which was based on 
the local density approximation, with the contact density for the homoge- 
neous system obtained by interpolating between the BEC, unitary, and BCS 
limits. These results provide direct experimental verification of the role of the 
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Fig. 6 Test of the adiabatic relation in a trapped gas of 40 K atoms, from Ref. |39) . The 
dimensionless contact C/kp as a function of l/kpa determined from the derivative of 
the energy E with respect to 1/a (solid dots) is compared to the measurements using 
photocmission spectrometry (open dots). The Fermi wavenumber kp is defined by the 
Fermi energy for the trapped system: Ep = h 2 kp/2m. 



contact in large-momentum and high-frequency tails for a many-body system 
of fermions with a large scattering length. 



5-4 Tests of the thermodynamic Tan relations 



The adiabatic relation in Eq. (|8| and the virial theorem in Eq. (11) relate 
different contributions to the total energy E to the contact C. The adiabatic 
relation expresses a derivative of (T + U) + V in terms of the contact. The 
virial theorem expresses the combination (T + U) — V in terms of the contact. 
Thus measurements of T + U, V, and C provide two nontrivial tests of the 
Tan relations for the thermodynamic properties of the system. 

The Jin group at JILA in Boulder has tested these Tan relations by mea- 
suring T+U and V for the same system of 40 K atoms for which they measured 
the contact C [35], as described in Section 5.3 They measured the external 
potential energy V by imaging the spatial distribution of the cloud of atoms, 
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which was trapped in a harmonic potential. They measured the combination 
T + U, which can be called the release energy, by turning off the trapping 
potential and observing the resulting expansion of the cloud. They measured 
T+U and V at values of l/A^a ranging from about —3 to about +0.3. For the 
contact C, they used their measurements from photoemission spectroscopy 
described in Section 5.3 They found good agreement between the two sides 
of the adiabatic relation in Eq. ^ as shown in Figure [6] They also found 
that the two sides of the virial theorem in Eq. (11) agreed to within the er- 
rors, which were roughly 1% of the Fermi energy. These results provide direct 
experimental verification of the role of the contact in the thermodynamics of 
a many-body system of fermions with a large scattering length. 



5.5 Contact density near unitarity 

The contact density for the homogeneous gas has been determined by the 
Salomon group at Ecole Normale Superieure in Paris [3D]. They used an 
imbalanced mixture of 6 Li atoms in the lowest two hyperfine spin states at 
a magnetic field near 834 G. They determined the equation of state for the 
homogeneous gas by measuring the number densities of the two spin states in 
a harmonic trapping potential integrated over the two tranverse dimensions. 



Their result for the numerical constant in Eq. (23) was C = 0.93(5). This is 



consistent with the value in Eq. (24) obtained from diffusion Monte Carlo 



calculations of the density-density correlator. 



6 Derivations of Universal Relations 

In this section, we give an overview of various derivations of the universal 
relations. We begin by describing briefly the novel methods used in the orig- 
inal derivations of the Tan relations. We then describe briefly various other 
approaches that have been used to rederive the Tan relations. Finally, we 
describe in more detail how universal relations can be derived using quan- 
tum field theory methods involving renormalization and the operator product 
expansion. 



6.1 Preliminaries 



The scattering amplitude for S-wave atom-atom scattering can be written as 
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~ kcotS (k)-ik' (49) 

where k is the relative wavenumber and So(k) is the S-wave phase shift. If 
the interactions have a finite range, the low-energy expansion of the phase 
shift can be expressed as a power series in k 2 : 

k cot 5 (k) = -l/a+\r s k 2 + (50) 

where r s is the effective range. The coefficient of (fc 2 )™ has dimensions 
(length) 2,1-1 . Generically, all these coefficients are comparable to the range 
raised to the appropriate power. By fine tuning the interactions, the scatter- 
ing length can be made much larger than the range. This fine-tuning leads to 
universal properties that depend on the interactions only through the scat- 
tering length. 

The universality for large scattering length reflects a well-behaved zero- 



range limit, in which all the terms on the right side of Eq. ( 50 ) go to zero 



except the leading term. The phase shift reduces in this limit to 

fccot<5 (fc) = -1/a. (51) 

The model in which the phase shift has this simple form up to arbitrarily 
large momentum is called the Zero-Range Model. The universal properties of 
a general model with large scattering length are realized in the Zero-Range 
Model in a particularly simple form, because a is the only length scale that 
arises from interactions. The price that must be paid is that the zero-range 
limit leads to divergences in some observables. It also leads to singularities 
in intermediate steps of the derivations of universal relations. An illustration 
is provided by the energy relation in Eq. The kinetic energy T and the 
interaction energy U are separately ultraviolet divergent, but the divergences 
cancel in the sum T + U. The singularities associated with the zero-range 
limit can be regularized by backing off from the zero-range limit or by some 
equivalent device. 

One way to represent the Zero- Range Model is in terms of the Schrodinger 
equation for noninteracting particles with peculiar boundary conditions. The 
Schrodinger wavefunction for Ni fermions in state 1 and N 2 fermions in state 
2 is a complex function <^(ri, . . . , rjv t ; r[, . . . , r' Nl ) that is totally antisym- 
metric in the first Ni positions and totally antisymmetric in the last N 2 
positions. The proper normalization of the wavefunction is 

Vxkl / d3ri ' ' ' S^™ 1 / dVl ' ' ' I d3r ' N2 l* (ri ' -' rjVi;ri ' = L 



(52) 



In the zero-range limit, this wavefunction diverges when the positions and 
r'j of any pair of fermions with different spins coincide. Its behavior when r% 
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and r[ are nearly equal is 

&(R+\r,r 2 , ...,Y Nl ;R-\r,r' 2 , ...,r' N2 ) 

— -> 4>{r) $(r 2 ,...,r Nl ;r' a ,...,r' N2 ;R), (53) 

where # is a smooth function of _R and </>(r) is the zero-energy scattering 
wavefunction for two particles interacting through a large scattering length: 

4>(r) = ---. (54) 
r a 

The Fourier transform of this wavefunction is 

ftk) = % - ^(fc). (55) 

The Schrodinger equation for interacting particles with a large scattering 
length reduces in the zero-range limit to the Schrodinger equation for non- 
interacting particles with the wavefunction constrained to satisfy the Bethe- 
Peierls boundary conditions in Eq. ( 53 ) . 



6.2 Tan's derivations 



Tan derived many of his universal relations by using generalized functions, 
or distributions, to deal with the singularities associated with the zero-range 
limit [H[3]. He introduced distributions A(k) and L(k) whose values at finite 
k are 



A{k) 
L(k) 



\k\ < oo, 
Ifcl < oo, 



and which have the following integrals over all k : 

d3k 1 ,/n n 

;A(k) = 0, 



(2tt) 3 k- 
d 3 k 1 



L(k) = 1. 



(56a) 
(56b) 



(57a) 
(57b) 



(2tt) 3 k 2 ' 

Using these properties, it is easy to verify that the Fourier transform of the 



zero-energy scattering wavefunction in Eq. ( 55 1 satisfies 

1 



d 3 fc 
(2^)3 



A(k) 



Ana 



L(k) 



<j)(k) = 0. 



(58) 
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This property allows the Bethe-Peierls boundary condition in Eq. ( 53 ) to be 
expressed as an equality rather than as a limit: 



A(k) 



1 

Attci 



-L(k) 



d 3 r e~ tk - 



X&(R+ir, r 2 , ...,r Nl ;R~\r, r' 2 , r' N ) = 0. 



(59) 



In addition to using A(k) and L(k) to impose the Bethe-Peierls boundary 
conditions, Tan used the distribution A{k) to regularize the interaction term 
in the Hamiltonian. He derived many of his universal relations by simple 
manipulations involving these distributions. They lead to an expression for 
the contact of the form 



2 ^ 



d 6 k 
(2tt)3 



L(k)k 2 n a (k). 



(60) 



By the properties of the distribution L(k) in Eqs. (56b) and (57b), the integral 
extracts the coefficient of 1/fc 4 in the high-momentum tail of n a (k). 



6.3 Other derivations 



6.3.1 Braaten, Kang, and Platter 



Braatcn and Platter [3] used quantum field theory methods to rcdcrivc all the 
Tan relations in Section [2] except for the relation in Eq. (10) for the sudden 



change in the scattering length. They used the formulation of the Zero- Range 
Model as a local quantum field theory, as described later in Section 6.4. 1| 



The singularities associated with the zero-range limit were regularized by 
imposing an ultraviolet cutoff on the momenta of virtual particles. The con- 
tact density was identified as the expectation value of a local operator. The 
universal relations were derived using renormalization, as described in Sec- 
tion 6.4.2[ and the operator product expansion, as described in Sections 6.4.4 
and f6A5l 

Braaten, Kang, and Platter [50] used quantum field theory methods to 
derive universal relations for the Resonance Model, in which the S-wave scat- 
tering phase shift is given by 

fccoa (fc) = - (\ + jf—^j ■ (6i) 

The Resonance Model is a 2-channel model, in which the states in the 2-atom 
sector consist of a point-like molecule as well as the usual 2-atom scattering 
states. It provides a natural model for a Feshbach resonance. The scatter- 
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ing length a = A — g 2 /v agrees with the expression in Eq. (37 1 if we set 
A = flb a , v = —mfi mo i(B — B Q )/h 2 , and g 2 — 1/i?*, where is given in 
Eq. (39 1. Braaten, Kang, and Platter found that in the various universal re- 
lations that correspond to the Tan relations, the contact density is replaced 
by various linear combinations of the expectation values of three local op- 
erators. If the expression for kcotSo(k) in Eq. (61 1 is well-approximated by 
— 1/a for all wavenumbers from up to the scale set by the system, the ex- 
pectation values of the three local operators must coincide in order for the 
universal relations to reduce to those of the Zero-Range Model. For a broad 
Feshbach resonance, which is defined by |fflbg| S> R*, this requires only that 
\a\ 3> ro and kr$ -C 1, where r$ is the range of interactions in the absence of 
the Feshbach resonance. For a narrow Feshbach resonance, which is defined 
by |abg| ^ R*> this requires also that kR* <C 1. 



6.3.2 Werner, Tarruell, and Castin 



Werner [57] derived the virial theorem in Eq. (11) independently. The right 



side of Eq. ( 11 ) was expressed not in terms of the contact C but in terms of 
the derivative of the energy that appears in the adiabatic relation in Eq. 

Werner, Tarruell, and Castin [35] rederived the tail of the momentum dis- 
tribution in Eq. ^ , the adiabatic relation, and the density-density correlator 
at short distances in Eq. Q . The singularities associated with the zero-range 
limit were regularized by using a lattice model in which the fermions oc- 
cupy the sites of a 3-dimensional cubic lattice whose spacing 6 approaches 
0. Werner, Tarruell and Castin also used a two-channel model to derive the 



universal relation for the number of closed-channel molecules in Eq. ( 38 ) . 

Werner and Castin [41] subsequently presented a much more thorough 
treatment of the universal relations for the Zero-Range Model and for the 
lattice model. In addition to the tail of the momentum distribution, the adi- 
abatic relation, and the density-density correlator at short distances, they 
rederived the energy relation in Eq. Q and the adiabatic sweep theorem, 



which is the generalization of Eq. ( 45 1 that allows for a time-dependent exter- 
nal potential. They showed that for a system in thermal equilibrium (which 
includes the ground state as a limiting case), the contact is an increasing 
function of 1/a: 

This inequality holds whether the derivative is evaluated at fixed entropy or 
at fixed temperature. The monotonic increase of the contact density with 
1/a is illustrated in Figs. [I] and [2] for the cases of a homogeneous gas at zero 
temperature that is balanced and strongly imbalanced, respectively. Because 
the inequality in Eq. ( [62] ) does not hold in general, it should not be regarded 
as a universal relation. 
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Werner and Castin also generalized the universal relations to various other 
systems [H] . They considered the effects of a nonzero range for the interaction 
potential. They derived the generalizations of the universal relations to two 
spatial dimensions, which had been considered previously by Tan [42] . They 
also derived universal relations for two types of fermions with unequal masses 
mi < m2, for identical bosons, and for mixtures of fermions and bosons 
with various masses. In some cases, including identical bosons and two types 
of fermions whose mass ratio m2/mi exceeds the critical value 13.7, the 
universal relations are complicated by the Efimov effect [H] . 



6.3.3 Zhang and Leggett 

Zhang and Leggett [H] rederived the adiabatic relation in Eq. d8| and the 



pressure relation in Eq. (12 1. They used a nonlocal quantum field theory with 
field operators i/ji(r) and ^> 2 (r) and with a short-range interaction potential 
U(r). They restricted their attention to a homogeneous many-body system 
in equilibrium. The equal-time 2-particle correlation function was expressed 
as a sum over eigenstates of the operator V ; 2( r 2)V'i( r i) : 

<^(ri)4(r 2 )^(r 2 )Vi(n)) = W (n,r 2 )* 0«(n,r a ), (63) 

i 

where v% is the average number of pairs with different spins in the eigenstate 
i and the eigenfunctions </>^(ri, r 2 ) are normalized accordingly in a large 
volume V. The Bethe-Peierls boundary conditions together with translation 
invariance imply that the limiting behavior of these eigenfunctions as the 
separation r = |n — r 2 | goes to zero is 

(j>®{R+\r,R-\r) — > C w e iP ^- R 0(r), (64) 

where the normalization constant CW and the momentum vector PW depend 
on the eigenstate i and 4>(r) is the zero-energy scattering wavefunction defined 



in Eq. ( 54 ) . The integral of the correlation function weighted by an arbitrary 



short-distance function s(\ri — r 2 |) therefore reduces to: 

1 f°° 

= —VC dr r 2 s{r)\(/)(r)\ 2 , (65) 



o 



47T 

where C is the contact density, which is given by 



C = 16tt 2 5> 1 |C«| 2 . 

i 



(66) 



Universal Relations for Fermions with Large Scattering Length 



31 



Zhang and Leggett emphasized that the thermodynamics is universal and 
is completely determined by the contact density C. Zhang and Leggett also 
used a two-channel model to derive a universal relation for the number of 



closed-channel molecules that can be reduced to Eq. ( 38 ) 29J . 

Zhang and Leggett derived a factorization formula for the interaction en- 
ergy density |29j that separates the dependence on a and thermodynamic 
variables, such as the temperature and number densities, from the depen- 
dence on the short-distance parameters that determine the shape of the in- 



teraction potential U (r). The interaction energy is the special case of Eq. (65 1 
in which the short-range function s(r) is the interaction potential U(r). The 
interaction energy density can be expressed as 

1 f°° 

U= — C dr r 2 U{r)Wr)\ 2 . (67) 
4tt Jo 

All the dependence on the thermodynamic variables and on the large scat- 
tering length a resides in the contact density C. 



6.3.4 Combescot, Alzetto, and Leyronas 

Combescot, Alzetto, and Leyronas E3] rederived the tail of the momentum 
distribution in Eq. ^ and the energy relation in Eq. |6]). They used the 
Schrodinger formalism in the coordinate representation to describe a sys- 
tem consisting of N% + N 2 fermions. The singularities associated with the 
zero-range limit were regularized by imposing a cutoff |r, — r'j\ > tq on the 
separations of the two types of fermions. They expressed the contact in terms 
of the Fourier transform of the function <P defined by the Bethe-Peierls bound- 



ary conditions in Eq. (53). By reverting to the coordinate representation, the 



contact density can be expressed as 



16tt 2 



(JVi - l)!(iV 2 - 1)! 



C ( R ) = 7Tr f^u TTT / d3 ^--- d 3 r Wl /dV 2 ... d 3 ~< 



r 



N 2 



1 1 2 

x \<P(r 2 ,...,r Nl ;r' 2l ...,r' N ;R)\ 



provided the wavefunction & is properly normalized as in Eq. ( [52] ) . This 
expression for the contact density was originally derived by Tan [S] - 

Combescot, Alzetto, and Leyronas generalized the tail of the momentum 
distribution and the energy relation to various other systems [43] . They gener- 
alized them to two spatial dimensions, which had been considered previously 
by Tan [32] • They generalized them to two types of fermions with different 
masses mi < rri2- They did not however consider the complications associ- 
ated with the Efimov effect when the mass ratio mijm\ exceeds the critical 
value 13.7 ED. 
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6.4 Quantum field theory derivations 

We proceed to explain how universal relations can be derived concisely using 
the methods of quantum held theory. 



6.4.1 Zero-Range Model 

A quantum held theory that describe atoms with two spin states must have 
fundamental quantum helds ip a (r), a — 1,2. The Hamiltonian operator for 
a local quantum held theory can be expressed as the integral over space of a 
Hamiltonian density operator: H = J RW. If the atoms are in an external 
potential V(r), the Hamiltonian density operator is the sum of a kinetic term 
T, an interaction term U, and an external potential term V: 

H = T + U + V. (69) 

In the quantum held theory formulation of the Zero-Range Model, the three 
terms in the Hamiltonian density operator arc 

T = £^VV4-V4 A) 0R), (70a) 



u 



^44**A A) m (70b) 

V = V(R)J2^lMR)- (70c) 

a 

For simplicity, we have set ft = 1. The superscripts (A) on the operators 



in Eqs. (70a I and (70b I indicate that their matrix elements are ultraviolet 
divergent and an ultraviolet cutoff is required to make them well dehned. For 
the ultraviolet cutoff, we impose an upper limit |fe| < A on the momenta of 



virtual particles. In the limit A — > oo, the Hamiltonian density in Eq. (69 1 



describes atoms with the phase shift given by Eq. (51 1 provided we take the 
bare coupling constant to be 

/ a\ 47ra . 

9 ^ A) = T^2aAT, (71) 

In Ref. 0, Braaten and Platter identified the operator that measures the 
contact density in the quantum held theory formulation of the Zero-Range 
Model. It is convenient to introduce the diatom held operator <P dehned by 

4>(R)=g (A)ip 2 4 A) (R). (72) 



This operator annihilates a pair of atoms at the point R. The superscript 
(A) on the operator V'iV'2 indicates that its matrix elements are ultraviolet 
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divergent. Their dependence on A is exactly compensated by the prefactor 
go(A), so <P is an ultraviolet finite operator. The contact density operator is 
[5]. This is just the interaction energy density operator multiplied by a 
constant that depends on the ultraviolet cutoff: <P^<P = mgoU. The contact 
is obtained by taking the expectation value of the contact density operator 
and integrating over space: 

C= [d 3 R (*t$(jj)). (73) 



6.4.2 Renormalization 



Several of the Tan relations follow very simply from the renormalization of 



the Zero- Range Model [5j . The renormalization condition in Eq. ( 71 1 implies 
that the bare coupling constant go(A) satisfies 



Its derivative with respect to a is 



(74) 



da 



9o(A) 



Ana 2 



9l{A). 



(75) 



The energy relation in Eq. ^ and the adiabatic relation in Eq. ([8| follow 
simply from these properties of the bare coupling constant: 

Energy relation. The kinetic and interaction terms T and U in the 



Hamiltonian density operator are given in Eqs. (70a) and (70b I. After 



inserting the expression in Eq. ( 74 ) for the bare coupling constant go into 



U, the sum of T and hi can be expressed as the sum of two finite operators: 



(76) 



By taking the expectation value of both sides of Eq. ( 76 1 , integrating over 



space, and using the expression for C in Eq. (73). we obtain the energy 
relation in Eq. 

Adiabatic relation. According to the Feynman-Hellman theorem, the 
rate of change in the energy due to a change in the scattering length can 
be expressed in the form 
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The Hamiltonian density % depends on a only through the factor of go 



in U. Using the derivative of the bare coupling constant in Eq. (75), we 
obtain the derivative of the Hamiltonian density: 

— = 1 <gt<j>. (78) 
da Airma 2 



By inserting this into Eq. ( 77 1 and using the expression for C in Eq. ( 73 1 



we obtain the adiabatic relation in Eq. n8h 



The virial theorem in Eq. (11) and the pressure relation in Eq. (12) can 
be derived by combining renormalization with dimensional analysis. For the 
purposes of dimensional analysis, we can regard h = 1 and m simply as 
conversion constants that allow any dimensionful quantity to be expressed as 
a length raised to an appropriate power. 

Virial theorem. For a system in a harmonic trapping potential, the only 
parameters that an energy eigenvalue can depend on are the scattering 
length a and the angular frequency w. The combinations with dimensions 
of length are a and (mw) -1//2 . Since an energy eigenvalue E has dimensions 
1/ (to length 2 ), the constraint of dimensional analysis can be expressed 
as the requirement that a differential operator that counts the factors of 
length gives —2 when acting on E = J d 3 R {%): 

a^- 2ui^—\ [d 3 R(W) = -2E. (79) 



da 8uj 

Using the Feynman-Hellman theorem, this equation can be written 



d 3 R ( — *— - 4(V) ) = -IE. (80) 



Airma 



Eq. (11) 



Using the expression for C in Eq. ( 73 1 , we obtain the virial theorem in 



Pressure relation. For a homogeneous system, the only variables that the 
free energy density J- = £ — TS can depend on are the scattering length 
a, the temperature T, and the number densities rii. The combinations 
with dimensions of length are a, (to/csT 1 ) -1 / 2 , and n i 1 ^ 3 . Since T has 
dimensions 1/ (to length 5 ), the constraint of dimensional analysis can be 
expressed as the requirement that a differential operator that counts the 
factors of length gives —5 when acting on T: 

-k-^lr-^k-^iL)^-^ (81) 

Using the adiabatic relation in Eq. mb, this can be written 
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—— C + 2TS - 3/Xim - 3/i 2 n 2 = -5J", (82) 
Anma 

where [ii is the chemical potential for the spin state i. The pressure relation 
in Eq. ( 12 I then follows from the thermodynamic identity T = —V+(iini + 



6.4.3 The operator product expansion 

Many universal relations are most concisely derived using the operator prod- 
uct expansion (OPE) of quantum held theory. The OPE was invented inde- 
pendently in 1969 by three giants of theoretical physics: Leo Kadanoff [33], 
Alexander Polyakov [35], and Ken Wilson [3B]. It is an expansion for the 
product of local operators at nearby points in terms of local operators at the 
same point: 

O a (R + |r) O b {R - |r) = /a.bW O c (R). (83) 

c 

The Wilson coefficients f% B (r) are ordinary functions of the separation vec- 
tor r. The local operators Oc{R) include some that can be obtained by 
Taylor-expanding the operators on the left side in powers of r, but they also 
include additional operators that take into account effects from quantum 
fluctuations. The Wilson coefficients for these operators are not necessarily 
analytic functions of the vector r, and they can even diverge as r — > 0. One 
particularly simple local operator is the unit operator /, whose expectation 



value in any state is 1. The sum over C in Eq. (83 1 can be extended to include 
all local operators if we allow Wilson coefficients that are 0. 

The OPE is a natural tool for generating universal relations, because it 
is an operator identity. By taking the expectation value of both sides of the 



OPE in Eq. (83) in some state of the system, we get an expression for the 
correlator of the operators Oa and Ob in terms of the expectation values of 
local operators in that state. Since this expansion holds for any state of the 
system, it is a universal relation. 

A local operator Oc (R) is assigned the scaling dimension dc if the corre- 
lation function of Oc and its hermitian conjugate at points separated by r 
decreases asymptotically as l/r 2dc at small r. The unit operator / is assigned 
scaling dimension 0. In a weakly-interacting theory, the scaling dimensions 
can be obtained simply by dimensional analysis. In a strongly-interacting 
theory, they can be significantly different. The difference between the scal- 
ing dimension and its value in the corresponding noninteracting theory is 
called the anomalous dimension. At very small r, the leading behavior of a 
nontrivial Wilson coefficient is determined by the scaling dimensions of the 
operators: 

/abW-^-^-*. (84) 
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In the OPE in Eq. ( 83 ) , the Wilson coefficients of higher dimension opera- 
tors go to more rapidly as r —> 0. A Wilson coefficient can be suppressed 
by a further power of r is there is an explicitly broken symmetry which, if 
exact, would require /^s( T ") t° vanish. The extra suppression factor is the 
dimensionless combination of r and the symmetry breaking parameter. 

The technical assumptions required to prove the OPE have been discussed 
by Wilson and Zimmerman [4 7) . The OPE can be expressed more precisely 



as an asymptotic expansion in the separation r = \r\. The OPE in Eq. (83 1 
can be organized into an expansion in powers of r by expanding the Wilson 
coefficients as Laurent series in r. The OPE is an asymptotic expansion if 
for any power p, there are only finitely many terms that go to zero more 



slowly than r p . The scaling behavior of the Wilson coefficients in Eq. (84 1 
guarantees that the OPE is an asymptotic expansion provided every local 
operator Oc(R) has a positive scaling dimension dc > and there are only 
finitely many local operators with scaling dimension dc < d for any positive 
number d. These conditions are satisfied in the Zero-Range Model, the Reso- 
nance Model, and other renormalizable local quantum field theories that are 
relevant to cold atoms. 

An illustration of the operator product expansion with anomalous dimen- 
sions is provided by the Ising Model in 2 dimensions. The exact solution by 
Lars Onsager in 1944 [48 implies that correlation functions in the contin- 
uum limit have scaling behavior with anomalous dimensions. For example, 
the leading term in the correlation function for two spin operators a as their 
separation r goes to has the form 

(o{R+\r)o(R-\r))^^ (85) 

where A is a constant that does not depend on the state of the system. This 
correlator is singular as r — > 0. The power law behavior suggests that the 
system becomes scale invariant at short distances. The fractional power of r 
indicates that the spin operator has an anomalous dimension. 



Kadanoff generalized the result for the correlator in Eq. ( 85 ) to an operator 
relation |44) : 

a(R + |r) a(R - §r) = ^1 + B\r\^e(R) + . . . , (86) 

where / is the identity operator, e(R) is the energy fluctuation operator, 
and B is another constant. The infinitely many terms that are not shown 



explicitly in Eq. (86 1 go to more rapidly than r 3 / 4 as r — >• 0. The short- 



distance tail of the correlator in Eq. ( 85 1 can be obtained simply by taking the 



expectation value of Eq. ( 86 ) . Kadanoff showed that the critical exponents 
of the Ising model, such as the exponent \ of \j\v\ in Eq. (85), could be 



deduced simply from the knowledge of which operator products nave singular 



Wilson coefficients. From the OPE in Eq. (86), we can deduce that the spin 
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operator has scaling dimension | and the energy fluctuation operator has 
scaling dimension 1. 



The OPE in Eq. (83) can be described more precisely as a short- distance 
operator product expansion, because it involves operators at the same time 
with small spatial separation. It can be generalized to a short-time operator 
product expansion, in which the operators also have a small separation in 
time: 

A {R+\r,T+\t) B {R-\r,T-\t)=Y,fA, B (r,t) O c (R,T). 

c 

(87) 

The Wilson coefficients f% B (r, t) are functions of the separation vector r and 
the time interval t. The short-time OPE is more subtle than the short-distance 
OPE because of the possibility that a correlator can have oscillatory behavior 
in t as t — > 47 . The possibility of oscillatory behavior is avoided in the 
Euclidean version of the quantum field theory that corresponds to analytic 
continuation of the time t to Euclidean time: t — > —it. Thus the short-time 



OPE in Eq. ( 87 ) can be expressed most rigorously as an asymptotic expansion 
in (r, —vr), where r is the Euclidean time separation obtained by the analytic 
continuation t — > —vr. In a Galilean-invariant theory, the appropriate scaling 
of (r, —vr) is such that r scales in the same way as \r\ 2 . 



6.4.4 Short-distance operator product expansion 

Universal relations for fermions with large scattering length can be derived 
from operator product expansions in the Zero-Range Model defined in Sec- 
tion |6.4.1| In this model, the scattering length a is the only length scale that 
arises from interactions. At distances much smaller than \a\, the model is scale 
invariant with nontrivial scaling dimensions. In the unitary limit a — > ±oo, 
the model is not only scale invariant at all distances but also conformally 
invariant |49j . If a is finite, we can regard 1/a as the symmetry- breaking 
parameter associated with the broken conformal symmetry. 

The fundamental field operators ipi and il>2 of the Zero-Range Model have 
the same scaling dimensions as in a noninteracting theory. However there 
are composite operators with anomalous scaling dimensions. The scaling di- 
mension of an operator Oc can be deduced from its propagator at large 
momentum k, which in a Galilean- invariant theory scales as k 2dc ~ 5 . Since 
the propagator for ipi or ?/>2 is simply (cj — k 2 /2m)~ 1 , these fields have scaling 
dimensions |. The scaling dimension of the number density operator ■f/'J.i/v 
is twice that of ip a , which is 3. If there were no interactions, the scaling 
dimension of the composite operator ipiif)2, or equivalently the diatom field 



operator <P defined in Eq. (72), would also be 3. However the propagator for 
<P is 
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fdt e iut Jd 3 r e- lk r ($(r,f) <? f (0,0)) 



-i47rm 



■1/a + \/—m(uj — k 2 /4m) 
Wc have dropped an additive constant that is independent of w and fc, which 



could be removed by renormalization. Since the propagator in Eq. ( 88 ) scales 
as fc -1 at large fc, <P has scaling dimension 2 and therefore anomalous dimen- 
sion — 1. The scaling dimension of the contact density operator is twice 
that of <P, which is 4. 

The short-distance OPE can be used to derive the Tan relation for the tail 
of the momentum distribution in Eq. (J2|. The momentum distribution n a (k) 
for atoms in the spin state a can be expressed as 

n a (k) = /d 3 i?|d 3 r e- ik - r {^l{R-\r) ^(R+^r)). (89) 

The behavior at large fc is dominated by the small-r region of the integral. 
We can therefore apply the OPE to the product of the operators ipl and ij) a . 
As shown in Ref. [5] , the leading terms in the OPE are 

i>l{R-\r) MR+\r)=^ a {R) 

-— *t*(Jl) + .... (90) 

We have written explicitly all terms whose Wilson coefficients go to zero more 
slowly than r 2 as r — > 0. The first two terms on the right side of the OPE 



in Eq. ( 90 ) can be obtained by multiplying the Taylor expansions of the two 
operators. The third term arises from quantum fluctuations involving pairs of 
atoms with small separations. That its Wilson coefficient is proportional to r 



can be predicted from the scaling dimensions of the operators using Eq. (84). 
The coefficient of r can be calculated using diagrammatic methods described 
in Ref. [30] . Note that this Wilson coefficient is not an analytic function of the 
vector r = (x, y, z) at r = 0, because it is proportional to r = \/ x 2 + y 2 + z 2 . 



The expectation value of the OPE in Eq. ( 90 ) can be expressed as 

(^t(fl_i r ) Vv(#+fr-)) = n a {R) + ir ■ V a (R) - ^ C ( R ) + ■ ■ ■ . 

(91) 

where T a is the momentum density of atoms in the state a. This form of the 
OPE was first written down by Tan [2]. 



Upon inserting the OPE in Eq. (90) into Eq. (89), the first two terms 



give a delta function in k and the gradient of such a delta function. They 
correspond to contributions to n a (k) that decrease at large fc faster than any 
power of fc. In the third term, the Fourier transform of the Wilson coefficient 
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at nonzero values of k can be obtained from the identity 

8tt 



d 3 r e~ lk - 



k 4 ' 



(92) 



which can be derived by differentiating the Fourier transform of a 1/r poten- 
tial. This term gives a power-law tail in the momentum distribution: 



n a {k) 



1 



d 3 R(&<P(R)). 



(93) 



Comparing with Eq. ^ , we obtain the expression in Eq. ( 73 ) for the contact 
C in the Zero-Range Model. This verifies that the contact density operator 
is &<P, where <l> is the diatom field operator defined in Eq. (72). 

The short-distance OPE can be used to derive the Tan relation in Eq. Q 
for the density-density correlator at short distances [5]. The OPE for the 
number density operators ipl^i an d VW2 includes a term whose Wilson co- 
efficient is singular as r — > 0: 



16tt 2 



ar J 



R) 



(94) 



All the other terms in the OPE are regular at r — 0. Taking the expectation 
value of both sides, we get the Tan relation in Eq. ([7|. 

The OPE in Eq. ( 94 ) can also be used to derive the universal relation for 



the static structure factor S\2(q) in Eq. (40). That structure factor can be 
expressed as 



S 12 (q) = J d3R f d3r e-^ r (^ 1 MR-hr)4MR+h r ))- 



(95) 



According to the universal relation in Eq. ( 40 ) , the high momentum tail 
has terms proportional to 1/q and 1/q 2 . They come from the singular terms 



proportional to 1/r 2 and 1 jr in the OPE in Eq. ( 94 ) 



6.4.5 Short-time operator product expansion 

Other universal relations can be derived using the short-time operator prod- 
uct expansion in Eq. (87). We will illustrate the use of the short-time OPE 
by showing how it can be used to derive the universal relations for radio- 
frequency (rf) spectroscopy that were presented in Section 



4.1 



The rf signal that causes a transition of an atom in spin state 2 into an 
atom in spin state 3 corresponds to the action of an operator - 03 - 02( r j £)• The 
inclusive rate -T(w) for the production of atoms in state 3 can be expressed 
as 



10 



r(«) = fi 2 Imijdt e^ +lt)t Jd 3 R Jd a 



(96) 



The symbol T in the matrix element implies that the product of operators is 
time ordered. For a large frequency w, the integrals are dominated by small 
Euclidean time intervals t = —it a nd b y small separations r. We can therefore 
apply the short-time OPE in Eq. (87 1 to the product of operators and 
ip\i^2- The Wilson coefficients of the leading one-body operator V4'02 an d the 
leading two-body operator were determined in Ref. |23| : 

Jdt e lut Jd 3 r ^ 3 (R+±r,t) ^^(R-^r, 0) 



- 4mr) 



-l 

12 



'13 )[ rt 12 



V- 



V- 



&<P(R) 



(97) 



Upon inserting the OPE in Eq. (97) into Eq. (96), the term gives a delta 
function in lu. This corresponds to a contribution that decreases at large oj 
faster than any power. The leading contribution to the high-frequency tail 
comes from the <pt<2> term: 



-l\2 



-Im 



-1 



j 1 - \J-m{u] + ie) 



'13 




R(&$(R)). 



(98) 



This reduces to the universal relation for the high-frequency tail in Eq. (33). 
Sum rules for the rf transition rate r(u>), such as those in Eqs. (31) and 



(32), can be derived by expressing the integral along the real to axis as a 



contour integral in the complex ui plane that wraps around the real axis. If 
we allow for a general weighting function f(ui), the sum rule becomes 



duj }(uj) r(u>) = n 2 



duj 
2^ 



f(u) dt e 



d 3 R ld 3 r 



x(T^ 3 (#+|r,i) 4MR-h,0)). 



(99) 



If that contour is deformed into a circle whose radius approaches inhnity, 
then u) has a large imaginary part along most of the contour. This justifies 
the use of the short-time OPE in Eq. (97). The sum rule can then be derived 



by evaluating the contour integral along the circle at infinity for each of 
the Wilson coefficients. For f(uS) = 1, the only nonzero contribution is from 
the - i/4' ! /'2 term and the sum rule reduces to Eq. (31). For f(uj) = u, the 



only nonzero contribution is from the <P^<P term and the sum rule reduces to 
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Eq. (32). The weighting function /(w) does not need to be a polynomial in 
oj. The Lorentzian function f(u>) — [(oj — loq) 2 + 7 2 ]" 1 gives a family of sum 
rules with two adjustable parameters that is less sensitive to range corrections 
than that in Eq. Q [25] . 

There are other universal relations that can be derived using the short- 
time OPE. One example is the high-frequency behavior of the structure factor 
Si2(uj,q), such as the universal relation in Eq. (41 1. It can be derived from 
the short-time OPE for the number density operators tp\ipi an d V4^2- 



Acknowledgements This research was supported in part by a joint grant from the Army 
Research Office and the Air Force Office of Scientific Research. I would like to acknowledge 
useful comments by Lucas Platter, Shina Tan, Edward Taylor, Felix Werner, and Shizhong 
Zhang. 



References 



1. E. Braaten and H. W. Hammer, Universality in Few-body Systems with Large Scat- 
tering Length, Phys. Rept. 428, 259 (2006), arXiv: cond-mat/0410417. 

2. S. Tan, Energetics of a strongly correlated Fermi gas, Annals of Physics 323, 2952 
(2008), arXiv:cond-mat/0505200. 

3. S. Tan, Large momentum part of a strongly correlated Fermi gas, Annals of Physics 
323, 2971 (2008), arXiv:cond-mat/0508320. 

4. S. Tan, Generalized virial theorem and pressure relation for a strongly correlated Fermi 
gas, Annals of Physics 323, 2987 (2008), arXiv : 0803 . 0841. 

5. E. Braaten and L. Platter, Exact relations for a strongly-interacting Fermi gas from the 
operator product expansion, Phys. Rev. Lett. 100, 205301 (2008), arXiv: 0803 . 1125. 

6. J.E. Thomas, J. Kinast, and A. Turlapov, Virial theorem and universality in a unitary 
Fermi gas, Phys. Rev. Lett. 95, 120402 (2005), arXiv : cond-mat/0503620. 

7. S. Tan, private communication. 

8. D. Blume and K.M. Daily, Universal relations for a trapped four-fermion system with 
arbitrary s-wave scattering length, Phys. Rev. A 80, 053626 (2009), arXiv : 0909 . 2701. 

9. S.Y. Chang, V.R. Pandharipande, J. Carlson, and K.E. Schmidt, Quantum 
Monte Carlo studies of superfluid Fermi gases, Phys. Rev. A 70, 043602 (2004), 
arXiv: physics/0404115. 

10. G.E. Astrakharchik, J. Boronat, J. Casulleras, and S. Giorgini, Equation of state of a 
Fermi gas in the BEC-BCS crossover: a quantum Monte Carlo study, Phys. Rev. Lett. 
93, 200404 (2004), arXiv : cond-mat/0406113. 

11. C. Lobo, I. Carusotto, S. Giorgini, A. Recati, and S. Stringari, Pair correla- 
tions of an expanding superfluid Fermi gas, Phys. Rev. Lett. 97, 100405 (2006), 
arXiv : cond-mat/0604282. 

12. D.S. Petrov, C. Salomon, and G.V. Shlyapnikov, Weakly bound dimers of fermionic 
atoms, Phys. Rev. Lett. 93, 090404 (2004), arXiv : cond-mat/0309010. 

13. R. Haussmann, M. Punk, and W. Zwerger, Spectral functions and rf response of ul- 
tracold fermionic atoms, Phys. Rev. A 80, 063612 (2009), arXiv : 0904 . 1333. 

14. Z. Yu, G.M. Bruun, and G. Baym, Short-range correlations and entropy in ultracold- 
atom Fermi gases, Phys. Rev. A 80, 023615 (2009), arXiv: 0905. 1836. 

15. F. Palestini, A. Perali, P. Pieri, and G.C. Strinati, Temperature and coupling depen- 
dence of the "contact" for an ultracold Fermi gas, arXiv : 1005 . 1158. 



42 



16. N. Prokof'cv and B. Svistunov, Fcrmi-polaron problem: Diagrammatic Monte Carlo 
method for divergent sign-alternating series, Phys. Rev. B 77, 020408(R) (2008), 
arXiv: 0707. 4259. 

17. M. Punk, P.T. Dumitrescu, and W. Zwerger, Polaron to molecule transition in a 
strongly imbalanced Fermi gas, Phys. Rev. A 80, 053605 (2009), arXiv: 0908. 1343. 

18. Z. Yu and G. Baym, Spin-correlation functions in ultracold paired atomic-fcrmion 
systems: Sum rules, self-consistent approximations, and mean fields, Phys. Rev. A 73, 
063601 (2006), arXiv: cond-mat/0510675. 

19. S. Zhang and A.J. Leggett, Sum-rule analysis of radio-frequency spectroscopy of ul- 
tracold Fermi gas, Phys. Rev. A 77, 033614 (2008). 

20. M. Punk and W. Zwerger, Theory of rf-Spectroscopy of strongly interacting fermions, 
Phys. Rev. Lett. 99, 170404 (2007), arXiv : 0707 . 0792. 

21. G. Baym, C.J. Pethick, Z. Yu, and M.W. Zwierlein, Coherence and clock shifts in 
ultracold Fermi gases with resonant interactions, Phys. Rev. Lett. 99, 190407 (2007), 
arXiv: 0707. 0859. 

22. M.W. Zwierlein, Z. Hadzibabic, S. Gupta, and W. Ketterle, Spectroscopic insensitivity 
to cold collisions in a two-state mixture of fermions, Phys. Rev. Lett. 91, 250404 (2003), 
arXiv : cond-mat/0306627. 

23. E. Braaten, D. Kang and L. Platter, Short-time operator product expansion for rf 
spectroscopy of a strongly interacting Fermi gas, Phys. Rev. Lett. 104, 223004 (2010), 
arXiv:1001.4518. 

24. W. Schneider, V.B. Shenoy, and M. Randcria, Theory of radio frequency spectroscopy 
of polarized Fermi gases, arXiv : 0903 . 3006. 

25. W. Schneider and M. Randcria, Universal short-distance structure of the single- 
particle spectral function of dilute Fermi gases, Phys. Rev. A 81, 021601(R) (2010), 
arXiv: 0910. 2693. 

26. P. Pieri, A. Perali, and G.C. Strinati, Enhanced paraconductivity-like fluctuations in 
the radiofrcquency spectra of ultracold Fermi atoms, Nature Physics 5, 736 (2009), 
arXiv: 081 1.0770. 

27. F. Werner, Virial theorems for trapped cold atoms, Phys. Rev. A 78, 025601 (2008), 
arXiv: 0803. 3277. 

28. F. Werner, L. Tarruell, and Y. Castin, Number of closed-channel molecules in the 
BEC-BCS crossover, Eur. Phys. J. B 68, 401 (2009), arXiv : 0807 . 0078. 

29. S. Zhang and A.J. Leggett, Universal properties of the ultracold Fermi gas, Phys. Rev. 
A 79, 023601 (2009), arXiv: 0809. 1892. 

30. E. Braaten, D. Kang, and L. Platter, Universal relations for a strongly- 
interacting Fermi gas near a Fcshbach resonance, Phys. Rev. A 78, 053606 (2008), 
arXiv: 0806. 2277. 

31. R. Combescot, S. Giorgini, and S. Stringari, Molecular signatures in the struc- 
ture factor of an interacting Fermi gas, Europhys. Lett. 75, 695 (2006), 
arXiv : cond-mat/0512048. 

32. H. Hu, X.-J. Liu, and P.D. Drummond, Static structure factor of a strongly correlated 
Fermi gas at large momenta, Europhys. Lett. 91, 20005 (2010), arXiv: 1003.3511. 

33. D.T. Son and E.G. Thompson, Short-distance and short-time structure of a unitary 
Fermi gas, Phys. Rev. A 81, 063634 (2010), arXiv: 1002 . 0922. 

34. E. Taylor and M. Randeria, Viscosity of strongly interacting quantum fluids: Spectral 
functions and sum rules, Phys. Rev. A 81, 053610 (2010), arXiv : 1002 . 0869. 

35. T. Enss, R. Haussmann, and W. Zwerger, Viscosity and scale invariance in the unitary 
Fermi gas, arXiv: 1008 . 0007. 

36. G.B. Partridge, K.E. Strecker, R.I. Kamar, M.W. Jack, and R.G. Hulet, Molecular 
probe of pairing in the BEC-BCS crossover, Phys. Rev. Lett. 95, 020404 (2005), 
arXiv : cond-mat/0505353. 

37. G. Veeravalli, E. Kuhnle, P. Dyke, and C.J. Vale, Bragg spectroscopy of a strongly 
interacting Fermi gas, Phys. Rev. Lett. 101, 250403 (2008), arXiv : 0809 . 2145. 



Universal Relations for Fermions with Large Scattering Length 



43 



38. H. Hu, E.D. Kuhnle, X.-J. Liu, P. Dyke, M. Mark, P.D. Drummond, P. Hannaford, 
and C.J. Vale, Universal behavior of pair correlations in a strongly interacting Fermi 
gas, Phys. Rev. Lett. 105, 070402 (2010), arXiv : 1001 . 3200. 

39. J.T. Stewart, J. P. Gaebler, T.E. Drake, and D.S. Jin, Verification of universal re- 
lations in a strongly interacting Fermi gas, Phys. Rev. Lett. 104, 235301 (2010), 
arXiv: 1002. 1987. 

40. N. Navon, S. Nascimbcnc, F. Chevy, and C. Salomon, The equation of state of 
a low-temperature Fermi gas with tunable interactions, Science 328, 729 (2010), 
arXiv: 1004. 1465. 

41. F. Werner and Y. Castin, Exact relations for quantum-mechanical few-body and 
many-body problems with short-range interactions in two and three dimensions, 
arXiv: 1001. 0774. 

42. S. Tan, S-wave contact interaction problem: a simple description, 
arXiv : cond-mat/0505615. 

43. R. Combcscot, F. Alzetto, and X. Leyronas, Particle distribution tail and related 
energy formula, Phys. Rev. A 79, 053640 (2009), arXiv: 0901 .4303. 

44. L.P. Kadanoff, Operator algebra and the determination of critical indices, Phys. Rev. 
Lett. 23, 1430 (1969). 

45. A. Polyakov, Properties of long and short range correlations in critical region, Sov. 
Phys. JETP 30, 151 (1970). 

46. K. G. Wilson, Non-Lagrangian models of current algebra, Phys. Rev. 179, 1499 (1969). 

47. K.G. Wilson and W. Zimmcrmann, Operator product expansions and composite field 
operators in the general framework of quantum field theory, Comm. Math. Phys. 24, 
87 (1972). 

48. L. Onsager, Crystal statistics. I. A two-dimensional model with an order-disorder tran- 
sition, Phys. Rev. 65, 117 (1944). 

49. D. T. Son and M. Wingate, General coordinate invariance and conformal invariancc 
in nonrclativistic physics: Unitary Fermi gas, Annals of Physics 321, 197 (2006), 
arXiv : cond-mat/0509786. 



